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Abstract 

In this paper, we consider the physics performance of a single far detector composed 
of a 100 kton next generation Liquid Argon Time Projection Chamber (LAr TPC) 
possibly located at shallow depth, coupled to the J-PARC neutrino beam facility with 
a realistic 1.66 MW operation of the Main Ring. The new far detector could be located 
in the region of Okinoshima islands (baseline L 658 km). Our emphasis is based 
on the measurement of the ^13 and 6cp parameters, possibly following indications for 
a non-vanishing ^13 in T2K, and relies on the opportunity offered by the LAr TPC 
to reconstruct the incoming neutrino energy with high precision compared to other 
large detector technologies. We mention other possible baselines like for example J- 
PARC-Kamioka (baseline L ~ 295 km), or J- PARC-Eastern Korean coast (baseline 
L ~ 1025 km). Such a detector would also further explore the existence of proton 
decays. 



*Now at IPHC, Universite Louis Pasteur, CNRS/IN2P3, Strasbourg, France. 
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1 Introduction 



The "Quest for the Origin of Matter Dominated Universe" is a long standing puzzle of our 
physical world (see e.g. [1]). Answers to this question might come from further exploration 
of 

• the Lepton Sector CP Violation by precise testing of the neutrino oscillation processes 

— measure precisely the 5c p and the mixing angle ^13; 

— examine matter effect in neutrino oscillation process and possibly conclude the 
mass hierarchy of neutrinos. 

• Proton Decay: 

— Search for p — z/ K"*" and p — > e 7r° in the range 10^^"^^^ years 

with assuming non-equilibrium environment in the evolution of universe. 

The primary motivation of the T2K experimenlQ [2] is to improve the sensitivity to the 
conversion phenomenon in the atmospheric regime by about an order of magnitude 
in the mixing angle sin^ IQn compared to the CHOOZ experimental limit |3j. The discovery 
of a non-vanishing ^13 angle would ascertain the 3x3 nature of the lepton flavor mixing 
matrix. The measurements in T2K might indicate that the ^13 angle is in a region where the 
simultaneous determination of neutrino mass hierarchy and the CP violating phase becomes 
possible. 

The J-PARC accelerator complex which includes the 180 MeV LINAC, the 3 GeV Rapid 
Cycling Synchrotron (RCS) and the 30-50 GeV Main Ring Synchrotron (MR) is planned to be 
commissioned in 2008. The J-PARC neutrino beam facility, under construction for the T2K 
experiment, is foreseen to begin operation in 2009. The final goal for the T2K experiment is 
to accumulate an integrated proton power on target of 0.75 MWx5 x 10^ seconds. Within 
a few years of run, critical information, which will guide the future direction of the neutrino 
physics, will be obtained based on the data corresponding to about 1^2 MW x 10'^ seconds 
integrated proton power on target (roughly corresponding to a 3(J discovery at sin^2^i3> 
0.05 and 0.03, respectively) [1]. 

It is well documented in the literature (see e.g. Refs. [5l[6l[7]) that the most challenging 
task for next generation long baseline experiments is to unfold the unknown oscillation 
parameters sin^ 2^13, bcp and mass hierarchy, sgn(Am3^), and from the measurement of the 
energy dependence of the oscillation signal to resolve the so-called problems of "correlations" 
and "degeneracies". 

The most important experimental aspects here are the beam profile (e.g. the ability to 
cover with sufficient statistics the 1^* maximum of the oscillation, the 1** minimum, and the 
2"'^ maximum), the visible energy resolution of the detector, with which the neutrino energy 

*We call T2K the approved Tokai-To-Kamioka experiment using Super-Kamiokande as far detector. New 
investigations or performance upgrades beyond this approved phase will presumably have a different name. 
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can be reconstructed, and the spectrum of the misidentified background (e.g. 7r° spectrum, 
etc.). 

One possible approach, advocated in the hterature [HI [9|, is to locate a very massive 
coarse Water Cerenkov detector at very long baselines {L > 2500 km) from a new multi- 
GeV neutrino source and to detect more than one oscillation maximum. In this kind of 
configuration, although the oscillation peaks are well separated in energy, an important 
issue is the ability to subtract the backgrounds, in particular that coming from misidentified 
neutral current events with leading vr^'s which introduce new sources of systematic errors. In 
order to alleviate this problem, a different approach assuming two similar far Water Cerenkov 
detectors, one located at Kamioka {L ~ 295 km) and the other in Korea {L ~ 1000 km) 
at the same OA2.5° from the J- PARC neutrino beam, was discussed in Ref. [lOj. The two 
detectors would see the same sub-GeV neutrino beam (in fact the same as in T2K) but 
the different baselines would allow to study E/L regions corresponding to the 1st and 2nd 
oscillation maxima. 

The approach of our paper is different. We consider the physics performance of a single 
100 kton next generation liquid Argon Time Projection Chamber (LAr TPC) located in the 
region of Okinoshima islands corresponding to a (relatively modest) baseline L ~ 658 km. 
A detector in this location will automatically see the J-PARC neutrino beam at an off-axis 
angle of ~ 0.8°. As we show below, this configuration allows to study both 1st and 2nd 
oscillation maximum peaks with good statistics. We rely on the realistic opportunity offered 
by the increase of intensity of the J-PARC MR and on the higher precision of the LAr TPC 
than other large detectors to separate the two peaks. In addition, the vr" background is 
expected to be highly suppressed thanks to the fine granularity of the readout, hence the 
main irreducible background will be the intrinsic z/g component of the beam. 

At this stage of our investigations, the emphasis is placed on the measurement of the 
6i3 and 6cp parameters, possibly following indications for a non-vanishing ^13 in T2K. We 
also compare the possibility to determine the sin^ 29i3-6cp parameters with a single detector 
configuration at other possible baselines like for example the T2K (Tokai-Kamioka, OA2.5°, 
L ~ 295 km) or the Tokai-Eastern Korean coast {L ~ 1025 km) with an off-axis OA1.0° pT| 
\T2\ . More complete physics performance studies will be presented elsewhere. 



2 A possible next step beyond T2K 

If a significant signal were to be observed at T2K, an immediate step forward to a next 
generation experiment aimed at the discovery of CP-violation in the lepton sector would be 
recommended with high priority. We have conducted a case study for the discovery of lepton 
sector CP violation based on the J-PARC MR power improvement scenario. Naturally, 
next generation far neutrino detectors for lepton sector CP violation discovery will be very 
massive and large. As a consequence, the same detector will give us the rare and important 
opportunity to discover proton decay. Thus, we also discuss the proton decay discovery 
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potential with a huge underground detector. 

Compared with the T2K experimental conditions, lepton sector CP violation discovery 
requires 

• an improved neutrino beam condition (intensity increase, broader energy spectrum, 
possibly re-optimization of the focusing optics, ...); 

• an improved far neutrino detector (by improvements we primarily mean increased signal 
reconstruction efficiency, better background separation, better energy resolution, 
and not only its volume). 

Detector improvements include 

• detector technology; 

• its volume; 

• its baseline and off-axis angle with respect to the neutrino source. 

As for the neutrino beam intensity improvement, a realistic first step power improvement 
scenario at J-PARC Accelerator Complex has been recently analyzed and proposed (LINAC 
energy is recovered to be 400 MeV, h=l operation at RCS and 1.92 seconds repetition cycle 
operation at MR). As for the detector technology, we assume a 100 kton Liquid Argon Time 
Projection Chamber for our case study. 

The effects of the CP phase 5c p appear either 

• as a difference between v and v behaviors (this method is sensitive to the CP-odd term 
which vanishes for 5cp = or 180°); 

• in the energy spectrum shape of the appearance oscillated charged current (CC) 
events (sensitive to all the non-vanishing 5cp values including 180°). 

Since antineutrino beam conditions are known to be more difficult than those for neutrinos 
(lower beam flux due to leading charge effect in proton collisions on target, small antineutri- 
nos cross-section at low energy, etc.), we concentrate on the possibility to precisely measure 
the V(. CC appearance energy spectrum shape with high resolution during a neutrino-only 
run. 



3 The liquid Argon Time Projection Chamber 

The liquid Argon Time Projection Chamber (LAr TPC) (See Ref. [13j and references therein) 
is a powerful detector for uniform and high accuracy imaging of massive active volumes. It 
is based on the fact that in highly pure Argon, ionization tracks can be drifted over distances 
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of the order of meters. Imaging is provided by position-segmented electrodes at the end of 
the drift path, continuously recording the signals induced. Tq is provided by the prompt 
scintillation light. 

In this paper, we assume a LAr TPC detector with a mass of the order of 100 kton, 
for example of a kind based on the GLACIER concept [14]. Unlike other liquid Argon 
TPCs operated or planned which rely on immersed wire chambers to readout the ionization 
signals, a double phase operation with charge extraction and amplification in the vapor phase 
is considered here in order to allow for very long drift paths and for improved signal-to-noise 
ratio. We assume that successful application of such novel methods will be an important 
milestone in the R&D for very large LAr TPC detectors in the range of 100 kton. At this 
stage, a ton-scale prototype based on this scheme has been developed and is under test |il5j . 
Challenges to realize liquid Argon TPCs with a scale relevant to this paper have been 
reviewed e.g. in Ref. jl6j . 

We refer to previous physics performance studies assuming such a detector configura- 
tion [121 Ej and recall that since Liquid Argon TPC has advantages on 

• good energy resolution/reconstruction, 

• good background suppression, 

• good signal efficiency 

it is suitable for a precision measurement of the neutrino energy spectrum to extract CP 
information. Thus we concentrate on the appearance energy spectrum shape measurement 
to extract leptonic CP phase information. 



4 The choice of far location: Okinoshima 

The J-PARC neutrino beam line was designed and constructed in order to allow an off-axis 
angle with respect to Super-Kamiokande (SK) between 2.5° and 3°. A beam setup yielding 
an OA2.5° at SK was chosen for the T2K experiment. In this configuration, the center of the 
T2K neutrino beam will go through underground beneath SK, and will automatically reach 
the Okinoshima island region with an off-axis angle ~ 0.8° and eventually the sea level east 
of the Korean shore with an off-axis angle ~ 1°. Larger off-axis angles are obtained moving 
inland Korea (either north, south or west). Figure [1] shows these baseline options using the 
same beam configuration as the T2K experiment. Parameters for different baselines and 
beam axes are summarized in Figure [H and Table [H 

The neutrino flavor oscillation probability including atmospheric, solar and interference 
terms, as well as matter effects, can expressed using the following equation [181 UHl [H] 

P{ue v^) ~ sin^ 2^13 ■ Ti + a ■ sin ^13 • (T2 + T3) + a" ■ T4. (1) 
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Okinoshima 



Japan^^^ ^^^^^^^^ ^ 




-1000km 
-1deg. Off-axis 



'658km 



295km 




'0.8deg. Off-axis jj 2.5deg. Off-axis 



Figure 1: Possible baseline candidates for 100 kton LAr TPC detector, Kamioka (short- 
baseline), Okinoshima (middle-baseline) and Korea (long-baseline). 



where, 



Ti = sin^0; 



23 ■ 



sin2[(l-A)- A] 
(T^Ap 

sin(AA) sin[(l - v4)A] 



T2 = sin bcp ■ sin IQyi ■ sin 26'23 ■ sin A 



T3 = cos bcp ■ sin 26^12 ■ sin 26*23 " cos A 



A {\-A) 
sin(AA) sin[(l - A)A] 



A {\-A) 

T4 = cos^^23-sin2 2^i2^^^^^. (2) 



where a = A = A = ^"^f^Pf . Amsi = - m?, Amsi = ml - ml is 

the mixing angle of the 1st and 3rd generations, while 6*12 is that for 1st and 2nd, and 6*23 is 
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Same beam configuration as T2K experiment 



Kamioka Okinoshima Korea eastern shore 



Baseline 
Off-axis angle 



295 km ~658 km ~1000 km 

2.5° ~0.8° ~1.0° 



Table 1: Parameters for baseline and axis with respect to beam of the candidates are 
summarized. 



that for 2nd and 3rd generations. 

The term which includes Ti is the "atmospheric term", those including T2 and T3 are 
"interference terms", and the one that includes T4 is the "solar term". The "interference 
terms" are sensitive to the 5c p phase and play therefore an important role to extract the 
CP phase. 

For definite calculations, we use the following parameters (we assume that most of these 
parameters will be precisely measured within the timescale of the one discussed here): 

• 6^23 is 7r/4. 

• is 0.572904 rad. 

• Am^^ is 8.2x10-5 eV^. 

• |Am2^| is 2.5x10^3 eV^. 

• Earth density for matter effects are 2.8 g/cm^. 

• normal hierarchy is assumed unless mentioned otherwise. 

Figure [2] shows the oscillation probability as a function of the E(GeV) / L(km) (neglecting 
for the moment matter effects). In the plot, the mixing angle sin^ 26'i3 was assumed to be 0.1. 
If the distance between source and detector is fixed, the curves can be easily translated to that 
for the expected neutrino energy spectrum of the oscillated events. As can be seen, if the 
neutrino energy spectrum of the oscillated events could be reconstructed with sufficiently 
good resolution in order to distinguish first and second maximum, useful information to 
extract the CP phase would be available even only with a neutrino run. 

We note that the following observables depend on leptonic CP phase: 

• position/height of the first oscillation maximum peak; 

• position/height of the second oscillation maximum peak. 

The position of the first oscillation minimum is unaffected by it (in this point all terms of 
the oscillation probability vanish except the solar term). To effectively experimentally study 
these observables, we point out that: 
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E(GeV) / L(km) 

Figure 2: Probability for oscillations as a function of the E(GeV)/L(km) for 

sin^ 2012, = 0.1 (for other oscillation parameters see text). Black curve shows the case for the 
6cp=0, red shows that for 6cp=90°, blue shows that for 6cp='270°, respectively. 



• the second oscillation maximum peak should end up at sufficiently high energy in order 
to be measurable; 

• the beam should have a sufficiently wide energy range to cover the 1st and 2nd maxi- 
mum; 

• the neutrino flux should be maximized in order to increase as much as possible the 
statistical signiflcance of the first and second maxima. 

In order to cover a wider energy range, we accordingly favor a detector location which is near 
on-axis. If one assumes that the second oscillation maximum has to be located at an energy 
larger than about 400 MeV, the oscillation baseline should be longer than about 600 km. In 
addition, in order to collect enough statistics, we choose a baseline which is not too much 
longer than above stated. 

Taking into account all of the above mentioned considerations, we privilege the Oki- 
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noshima region: placing a detector in an appropriate location on the island will probe 
neutrino oscillations at a baseline of ~ 658 km away from the source at an off-axis angle of 
- 0.8°. 



5 Neutrino Flux and expected Event Rates 

We assume realistic parameters of the J-PARC beam after all accelerator complex up- 
grades described in the KEK road-map are accomplished, as follows |1]: 

• The average beam power is reaching 1.66 MW; 

• A total of 3.45 x 10^^ POT is delivered on target per year; 

• The optimal kinetic energy of the incident protons is 30 GeV. 

To remain conservative, we focus on an analysis which uses a neutrino run only during 
five years under the best J-PARC beam assumption. An antineutrino beam (opposite horn 
polarity) might be considered in a second stage in order to cross-check the results obtained 
with the neutrino run (in particular see Section [8] on the mass hierarchy problem). The 
parameters of the Okinoshima location and the assumed beam are therefore as follows: 

• Distance from J-PARC is 658 km; 

• The axis of the beam is off by 0.76°; 

• five years operation with horn setting to neutrino run. 

The expected neutrino flux calculated under these assumptions is shown in Figure [3] 
where the curves correspond to one year run (3.45x10^^ POT). The black, red, green, blue 
lines show z/^, i/^, z/g, fluxes, respectively. 

The interacting neutrino cross section on Argon was computed using the NUANCE pro- 
gramme [20]. We use the foUowings parameters for the cross section calculation: 

• Number of protons is 18, and that of neutrons is 22. 

• Medium density is 1.4 g/cm^. 

• Nucleus Fermi Gas model with binding energy of 30 MeV and Fermi momentum of 
242 MeV/c. 

^The exact location of the potential experiment in Okinoshima has not yet been investigated, but we note 
that the island has several hills. 
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Figure 3: Calculated neutrino flux under the 

assumptions (3.45xl02i POT). Black shows pigure 4: Calculated neutrino- Argon cross 

z/^, red shows z/^, green shows u„ blue shows section, unit by pb. Nuclear effects are taken 

Ue fluxes, respectively. into account. 

Taking these into account, we obtain the neutrino cross section shown in Figure HI Table [2] 
shows that total number of charged current (CC) events at Okinoshima for the null oscillation 
case and the number of CC Ue events from i/^ v^. oscillations for three different mixing 
angles, and in various bcp scenarios, normalized to flve years neutrino run at 1.66 MW and 
100 kton fiducial mass. 

The number of expected events is already a good indicator of the signal, but for this 
analysis we use binned likelihood fitting for the energy spectrum, therefore all information 
including total number of events is taken into account, as described in the next section. 



6 Expected Energy Spectrum of z/g CC events 

Images taken with a liquid Argon TPC are comparable with pictures from bubble cham- 
bers. As it is the case in bubble chambers, events can be analyzed by reconstructing 3D-tracks 
and particle types for each track in the event image, with a lower energy threshold of few 
MeV for electrons and few tens of MeV for protons. The particle type can be determined 
from measuring the energy loss along the track [dE/dx) or from topology (i.e. observing 
the decay products). Additionally, the electronic readout allows to consider the volume as a 
calorimeter adding up all the collected ionization charge. The calorimetric performance can 
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Events for 100 kton at Okinoshima normalized to 5 


years 


at 1.66 MW beam power 










Beam components (null oscillation) 


82000 750 


1460 


35 


i^f. - 


^ Vf. oscillations 






6cp = 


0° 90° 


180° 


270° 


sin^2^i3 = 0.1 


2867 2062 


2659 


3464 


sin^ 2^13 = 0.05 


1489 1119 


1342 


1908 


sin^ 2^13 = 0.03 


942 506 


829 


1266 



Table 2: Number of CC events at Okinoshima: beam components for null oscillation, and 
oscillated events in various 8cp scenarios (normalized to five years neutrino run) 



be excellent, as we will show in Section [TOl depending on event energy and topology. 

In order to understand the effect of resolution on physics performance, we show in this 
section the z/g charged current (CC) event energy spectra using different simple models. One 
is the perfect resolution reference, and others are 100 MeV/200 MeV Gaussian 

resolution cases for more realistic cases. This "resolution" should include neutrino interac- 
tion effects as well as detector resolution. Possible smearing or backgrounds affecting the 
measurement of the neutrino energy spectrum are listed below subdivided into 4 classes: 

• Neutrino interaction: 

— Fermi motion and nuclear binding energy, 

— Nuclear interactions of final state particles within the hit nucleus (FSI), 

— Vertex nuclear remnant effects (e.g. nuclear break-up signal), 

— Neutral Current (NC) 7r° event shape including vertex activity. 

• Detector medium: 

— Ionization processes, 

— Scintillation processes, 

— Correlation of between amount of charge and light, 

— Charge and light quenching, 

— Hadron transport in Argon and secondary interactions, 

— Charge diffusion and attenuation due to impurity attachment. 

• Readout system including electronics system: 

— signal amplification or lack thereof. 
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— signal-to-noise ratio, 

— signal shaping and feature extraction. 

• Reconstruction: 

— Pattern recognition 

— Background processes (NC 7r°, CC, ...) and their event shape 

— Particle identification efficiency and purity 

Very few real events in liquid Argon TPCs have so far been accumulated to allow a full 
understanding of these complex effects and their interplay. The only sample comes from a 
small 50 It chamber developed by the ICARUS Collab. and exposed to the CERN WANF 
high energy neutrino beam which collected less than 100 quasi-elastic events [21j. 

It is clear that the energy resolution will depend on several detector parameters, including 
the readout pitch, the readout method chosen and on the resulting signal-over-noise ratio 
ultimately affecting the reconstruction of the events. We therefore stress that significantly 
improved experimental studies with prototypes exposed to neutrino beams of the relevant 
energies and sufficient statistics are mandatory to assess and understand these effects. In 
the meantime, we give in Section [10] preliminary estimates for potentially achievable energy 
resolutions in a liquid Argon medium. The results are based on full GEANT3 \2f2\ simulation 
of the energy deposited by final state particles in the detector volume, however do not include 
all possible contributing effects. 

Figure [5] shows the energy spectra of electron neutrino at the cases of 5c p equal 0°, 90°, 
180°, 270°, respectively. Shaded region is common for all plots and it shows the background 
from beam z/g. Here perfect resolution is assumed for reference to later cases. 

If we smear the energy spectra shown in Figure [5] with Gaussian of sigma equals to 100 or 
200 MeV independent from original neutrino energy, we obtain spectra shown in Figures [6] 
and [71 As seen easily, an energy resolution below 100 MeV is crucial since the robustness of 
the neutrino oscillation is directly determined by the visible second oscillation peak around 
400 MeV in the energy spectrum. In 200 MeV resolution case, the second peak is hidden by 
the smearing of the 1st oscillation maximum peak. 

7 Oscillation parameters measurement from energy 
spectrum analyses 

Assuming all others were measured very precisely, there are only two free parameters 
sin^26'i3 and 5cp to be fitted using the energy spectra. As shown in Figure [21 the value of 
5cp varies the energy spectrum, especially the first and the second oscillation peaks (heights 
and positions), therefore comparison of the peaks determine the value 5cPi while the value 
of sin^ 26'i3 changes number of events predominantly. 

To fit the free parameters, a binned likelihood method is used. The Poisson bin-by-bin 
probability of the observed data from the expected events (for an assumed pair of values 
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Figure 5: Energy spectra at sin^ 2^i3=0.03 
case, but 5cp = 0° (top-left), 90° (right-top), 
180° (left-bottom), 270° (right-bottom) cases. 




Neutrino Energy (GeV) 



Neutrino Energy (GeV) 



Figure 6: Energy spectra assuming Gaussian 
100 MeV smearing. 6cp = 0° (top-left), 90° 
(right-top), 180° (left-bottom), 270° (right- 
bottom) cases. 



(sin^ 26^13, Sep)) is calculated. The fit procedure is validated by testing the result on a 
pseudo-experiment. Figure [8] shows one typical pseudo-experiment for the perfect resolution 
case (with 6cp equals to 0, and sin^ 26^13 equals to 0.1). Best fit gives reasonable result. 

At this stage, we only take statistical uncertainty into account for the fitting, thus other 
systematic uncertainties like far/near ratio, beam z/g shape, energy scale, and so forth are 
not considered. Also, oscillated signal and beam z/g are only accounted in the fit, i.e. other 
background like neutral current vr*^ or beam z/g contamination is assumed to be negligible 
compared to beam i/^ contamination. 

As a reference, we first extract allowed regions in the perfect resolution case (See Fig- 
ure [9l). Twelve allowed regions are overlaid for twelve true values, sin^ 26'i3=0.1, 0.05, 0.02, 
and 5cp=0°, 90°, 180°, 270°, respectively. The 5cp sensitivity is 20~30° depending on the 
true 6c p value. 

Allowed regions are then extracted for 100 and 200 MeV resolution cases. The used 
energy spectra are same as Figures [6] and [71 Results are shown as Figures [TO] and [HI 

One obvious but important issue to be pointed out is the robustness of the fitting. The 
fit procedure shows that results could also be extracted with the 200 MeV resolution: this 
result is as expected statistically; however, we stress that in this case there is no second 
oscillation maximum peak visible in Figure [71 Hence, we think it is mandatory to keep an 
energy resolution less than 100 MeV as goal for the credibility of this experiment. 
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Neutrino Energy (GeV) 



Neutrino Energy (GeV) 




Figure 7: Energy spectra assuming Gaussian 200 MeV smearing. 5cp = 0° (top-left), 90° 
(right-top), 180° (left-bottom), 270° (right-bottom) cases. 



8 Investigation of mass hierarchy 

The influence of matter on neutrino oscillations was first considered by Wolfenstein |23j . 
As is well-known (although never directly experimentally verified), oscillation probabilities 
get modified under these conditions. Matter effects are sensitive to the neutrino mass or- 
dering and different for neutrinos and antineutrinos. As mentioned earlier, we consider in 
this paper the possibility of a neutrino-only run. Hence, we briefly address in this section 
the question of normal hierarchy (NH) versus inverted hierarchy (IH). Since we are focusing 
on the potential discovery of CP-violation in the leptonic sector, our discussion is geared 
towards possible ambiguities that would arise if the mass hierarchy was unknown. 

Indeed, Figure [12] illustrates the results of a fit of a pseudo-data with NH by both NH 
(black) and IH (red) hypotheses, assuming only the neutrino run. The best fit likelihood 
value with one assumed hierarchy is used to calculate the likelihood variation AL for both 
hierarchies. One could claim that CP-violation is discovered if the experimental results of a 
given experiment exclude the 5cp phase to be either or 180°. Hence, the danger is that the 
lack of knowledge of the mass hierarchy (or rather the "wrong" choice in hypothesis when 
selecting the hierarchy in the fit of the data), gives a result for 5cp consistent with either or 
180°. On the other hand, if fits of a given experiment with both assumed hierarchies provide 
neither nor 180°, one can declare "discovery" although the mass hierarchy could not be 
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Figure 8: Typical one pseudo-experiment to show the vahdity of the fitting. Histogram 
shows the best-fit (best oscillation shape -|- background), and the crosses show the pseudo- 
data. 



determined. Alternatively, if it turned out in the actual experiment that one of results of 
the fit is consistent with 6cp = or 180°, one would still have the possibility to consider an 
anti-neutrino run in order to solve this ambiguity. 

Our results indicate that CP-violation would be unambiguously discovered under our 
assumption for true values of (5cp's in the region of 90° and 270°. Alternatively, in case 
of true values of 6cp^s near 0° or 180° an antineutrino flux would help untangle the two 
solutions if the neutrino mass hierarchy was unknown. 



9 Expected sensitivity to sin^ 2^13 

Although this paper has focused on the possibility to measure the sin^ 26*13 and 6cp 
oscillation parameters, it is instructive to estimate the sensitivity of the potential setup in 
the case of negative result from T2K. The corresponding sensitivity to discover ^13 in the 
true (sin^ 26*13, (5cp) plane at 90% C.L. and 3a is shown in Figure fT3l where we assumed 5 
years of neutrino run comparing the 100 kton LAr at Okinoshima with increased statistics 
in SK at Kamioka and a potential 20 kton LAr at Kamioka. 

In order to discover a non- vanishing sin^ 26^13, the hypothesis sin^ 2^13 = must be 
excluded at the given C.L. As input, a true non-vanishing value of sin^ 26*13 is chosen in the 
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Figure 9: Allowed regions in the perfect res- 
olution case. Twelve allowed regions are over- 
laid for twelve true values, sin^ 2^i3=0.1, 0.05, 
0.02, and (5cp=0°, 90°, 180°, 270°, respec- 
tively. 



Figure 10: Allowed regions in the 100 MeV 
resolution (Gaussian sigma) case. Twelve al- 
lowed regions are overlaid for twelve true val- 
ues, sin^ 2^13=0.1, 0.05, 0.02, and (5cp=0°, 90°, 
180°, 270°, respectively. 



simulation and a fit with sin^ 2^13 = is performed, yielding the "discovery" potential. This 
procedure is repeated for every point in the (sin^ 26^13, 5cp) plane. 

At the 3(7 CL. the sensitivity of the T2K experiment is 0.02 [2J. Continuing to collect 
data at SK with an improved J-PARC neutrino beam for another 5 years would improve 
the T2K sensitivity by a factor ~ 2. In comparison, our simulations indicate that a 20 kton 
LAr TPC is expected to perform better than this although the masses are comparable, since 
the signal efficiency is higher than SK and the NC background is assumed to be negligible 
contrary to SK. Even better, a 100 kton LAr TPC at Okinoshima would further improve 
the sensitivity by about a factor six compared to SK at Kamioka (or a factor 10 compared 
to T2K), thanks to its bigger mass, increased cross-section at higher energies and reduced 
off-axis angle, although the neutrino fluxes at the same off-axis angle would be reduced by 
a factor ^ 5 because of the longer baseline. 
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Figure 11: Allowed regions in the 200 MeV resolution (Gaussian sigma) case. Twelve 

allowed regions are overlaid for twelve true values, sin^ 26'i3=0.1, 0.05, 0.02, and (5cp=0°, 
90°, 180°, 270°, respectively. 

10 Neutrino energy resolution in a LAr medium 

As mentioned previously, the neutrino energy resolution expected in LAr medium depends on 
several parameters and their interplay that ultimately need to be measured experimentally. 

In principle, energy and momentum conservation allow to estimate the incoming neutrino 
energy in the detector via a precise measurement of decay products, with the exception of 
the smearing introduced by Fermi motion and other nuclear effects (nuclear potential, re- 
scattering, absorption, etc.) for interactions on bound nucleons. In this section, we try to 
estimate the smearing introduced by these effects with the help of exclusive neutrino event 
final states distributed with the Okinoshima flux and generated with the GENIE MC |24] 
subsequently fully simulated with GEANT3. 

Nuclear effects in neutrino interactions can be roughly divided into those of the nuclear 
potential and those due to reinteractions of decay products. Bound nucleons and other 
hadrons in nuclei are subject to a nuclear potential. The Fermi energy (or momentum) must 
be calculated from the bottom of this nuclear potential well, and the removal of a nucleon 
from any stable nucleus is always an endothermic reaction. When hadrons are produced in 
the nucleus, some energy is spent to take it out of this well: for a nucleon, the minimum 
energy is given by the nucleon separation energy (around 8 MeV), and corresponds to a 
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Figure 12: Mass hierarchy investigation with neutrino run only. If fits with both hierarchy 
hypotheses provide neither nor 180°, one can declare discovery of CP violation in the 
leptonic sector. If any of the fits results in a of or 180°, then an anti-neutrino run 
could be envisaged. 
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Figure 13: Sensitivity to sin^(2^i3) for 5 years of neutrino run comparing 100 kton LAr at 
Okinoshima with SK at Kamioka and 20 kton LAr at Kamioka. 



nucleon at the Fermi surface. In this case, the daughter nucleus is left on its ground state. 
More deeply bound nucleons, leaving a hole in the Fermi sea, correspond to an excitation 
energy of the daughter nucleus, and an additional loss of energy of the final state products. 
This energy is then spent in evaporation and/or gamma deexcitation. Thus, the energy of 
the final state products is expected to be always slightly smaller than for interactions on free 
nucleons, and spread over a range of about 40 MeV. Correspondingly, the Fermi momentum 
is transferred to the decay products and compensated by the recoil of the daughter nucleus. 
Additional momentum distortions come from the curvature of particle trajectories in the 
nuclear potential. 

Reinteractions in the nuclear medium also play an important role. Interaction products 
can lose part of their energy in collisions, or even be absorbed in the same nucleus where they 
have been created. This is particularly true for pions, that have an important absorption 
cross section on nucleon pairs, while kaons have smaller interaction probability. 

Once final state products have exited the nucleus, they will propagate in the medium 
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Figure 14: Full GEANT simulations of the reconstructed neutrino energy spectra from 
detailed event simulation (see text) for neutrino oscillations with Sep = 0°: perfect recon- 
struction (top- left), LAr using deposited energy (right-top), ibid but with charged pion mass 
correction (left-bottom) using final state lepton only (like in WC) (right-bottom) cases. 



with the possibility that further interactions occur. 

In the case of the liquid Argon TPC, the medium is fully homogeneous and the detector 
is fully active. All deposited energy in the medium (above a certain threshold) will be 
eventually collected. Several methods can be adopted to reconstruct the neutrino energy 
and we list in the following three: 

1. Final state lepton: this method, traditionally used in large Water Cerenkov detectors 
like Kamiokande, 1MB or SK, relies on the precise measurement of the energy and 
direction of the outgoing lepton and kinematically constrains the incoming neutrino 
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energy by assuming a quasi-elastic configuration: 

^ ME, - m\l2 

^ M — Ei + Pi cos 9i 

where Ei, pi and cos 6£ are the energy, momentum and scattering angle of the outgoing 
lepton. This method is sensitive to the final state configuration and to Fermi motion 
(up to ~ 240 MeV/c) which randomizes the direction of the outgoing lepton. 

2. Momentum conservation: neglecting the incoming neutrino mass, one obtains 

h 

where pi is the momentum of the outgoing lepton, ph is the momentum of the outgoing 
hadron h and Pp is the Fermi motion of the hit nucleon (up to ~ 240 MeV/c). This 
method is also sensitive to Fermi motion since the recoiling remnant hit nucleus (with 
momentum —Pp) is not measured 0. 

3. Energy conservation: using a calorimetric approach, one can sum the deposited 
energies of all outgoing particles. In a tracking-calorimeter this is obtained by summing 
the dE/dx measurements along each ionizing track to obtain the associated kinetic 
energies T = J[dE/dx)dx. One should identify final state particles in order to take 
into account their rest masses. This method is sensitive to the nuclear binding energy 
of the decay products and is intrinsically more precise than the above method relying 
on momenta for the energies considered here. Mathematically one can write this result 
as: 

E, = Etot- M = Et + ^Eh- M = Et^^{Th + mh)- M 

h h 

= E, + T^ + Y,{T^±+m^±) + Y.E^ + ... (3) 



All methods are sensitive to potential re-interactions of the outgoing hadrons within the 
nucleus. In order to estimate the potential of a fully homogeneous and sensitive medium 
like the liquid Argon TPC, we performed full simulations of neutrino interactions: the event 
4-vector are generated with the GENIE MC and final state particles (after nuclear reinter- 
action) were propagated through the liquid Argon medium with GEANT3. The geometry of 
the setup is an infinite LAr box and at this stage we have neglected charge quenching (i.e. 

■''We point out that if one uses the knowledge on the direction of the incoming neutrino, the total (missing) 
transverse momentum can be kinematically constrained to zero, thereby providing a handle compensate for 
the transverse component of the Fermi motion. The longitudinal component of the momentum, which is not 
negligible at low neutrino energies, remains however undetermined. 
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we assume a linear response to dE/dx). We also assume a 100% efficiency to identify final 
state charged pions. 

The achievable incoming neutrino energy resolution in the liquid Argon medium has 
been estimated using v^. CC events distributed with the flux expected at the Okinoshima 
location. The results, using the energy conservation method, are shown in Figure [TH 



11 Other baseline options 

Although Okinoshima looks a very good candidate as mentioned so far, the favorable 
geography would allow, in principle, for a few baseline candidates for the detector locations; 
Kamioka for a relatively "short" baseline, Okinoshima for a medium baseline and Korea for 
the longest baseline. It is worthwhile to compare the physics potential at Okinoshima with 
other scenarios with the same strategy and same assumed beam conditions, i.e. five years 
neutrino run at 1.66 MW. 

Figure [15] illustrates the analyzed allowed regions using each option. A perfect resolu- 
tion was assumed here to directly compare the physics potential of the various sites. As 
seen, sensitivity of 5c p is similar between Okinoshima and Korea, but that in Kamioka is 
much worse than others. This indicates the analysis strategy using only neutrino run is not 
suitable for Kamioka case, as was expected since the OA2.5° available at Kamioka does not 
allow to cover 1st and 2nd oscillation maxima peaks. As for the sensitivity of the sin^2^i3, 
Okinoshima is better than Korea, due to the increased statistics at the shorter baseline. 
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Figure 15: Allowed regions if the location of the 100 kton LAr TPC were at Kamioka, 
Okinoshima or Korea, assuming the T2K beam optics is kept. Perfect resolution is assumed 
to illustrate the physics potential from each site selection. 
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12 Proton Decay Discovery Potential 



Grand Unification of tlie strong, weak and electromagnetic interactions into a single unified 
gauge group is an extremely appealing idea [25l [26] which has been vigorously pursued 
theoretically and experimentally for many years. The detection of proton or bound-neutron 
decays would represent its most direct experimental evidence. The physics potentialities of 
very large underground Liquid Argon TPC was recently carried out with detailed simulation 
of signal efficiency and background sources, including atmospheric neutrinos and cosmogenic 
backgrounds [27]. It was found that a liquid Argon TPC, offering good granularity and energy 
resolution, low particle detection threshold, and excellent background discrimination, should 
yield very good signal over background ratios in many possible decay modes, allowing to 
reach partial lifetime sensitivities in the range of 10'^^ — 10'^^ years with exposures up to 
1000 ktonxyear, often in quasi-background-free conditions optimal for discoveries at the 
few events level, corresponding to atmospheric neutrino background rejections of the order 
of 10^. Multi-prong decay modes like e.g. p —>■ fi~7T~^K~^ or p —>■ e^TT^Ti~ and channels 
involving kaons like e.g. p K^u, p — > e^K^ and p ^^K^ are particularly suitable, since 
liquid Argon imaging provides typically an order of magnitude improvement in efficiencies 
for similar or better background conditions compared to Water Cerenkov detectors. Up to 
a factor 2 improvement in efficiency is expected for modes like p —>■ e"'"7 and p —>■ fi'^'j 
thanks to the clean photon identification and separation from vr*^. Channels like p e^7r° or 
p /x'''7r°, dominated by intrinsic nuclear effects, yield similar efficiencies and backgrounds 
as in Water Cerenkov detectors. Thanks to the self-shielding and 3D-imaging properties of 
the liquid Argon TPC, the result remains valid even at shallow depths where cosmogenic 
background sources are important. In conclusion, a LAr TPC would not necessarily require 
very deep underground laboratories even for high sensitivity proton decay searches. 

Once again, we stress the importance of an experimental verification of the physics po- 
tentialities to detect, reconstruct and classify events in the relevant GeV energy range. This 
experimental verification will require the collection of neutrino event samples with high statis- 
tics, accessible e.g. with a detector located at near sites of long baseline artificial neutrino 
beams. 



13 Conclusions 

In this paper, we have reported our case study for a new future long baseline neutrino 
and nucleon decay experiment with a 100 kton Liquid Argon TPC located in the region of 
Okinoshima using the J-PARC Neutrino Facility. 

Our study assumes a realistic upgrade of the J-PARC Main Ring operation yielding an 
average 1.66 MW beam power for neutrinos coupled to the new far detector at Okinoshima, 
which would then see the beam at an off-axis angle of ~ 1° and at a baseline of ~ 658 km. 
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In this first discussion, we focused on tlie possibility to measure tlie 6*13 and 5c p pa- 
rameters. Our strategy is based on tlie precise measurement of the energy spectrum of 
the oscillated events, and in particular on the comparison of the features of 1st maximum 
oscillation peak with those of the 2nd maximum oscillation peak. In order to efficiently 
determine and study those features, we rely on the excellent neutrino event reconstruction 
and incoming neutrino energy resolution, as is expected in the case of the fully-sensitive and 
fine charge imaging capabilities of the liquid Argon TPC 

Beyond this case study, we intend to perform more detailed simulations of the perfor- 
mance of the experiment including for example different sources of systematic errors and 
sources of backgrounds which have been up-to-now neglected. 

The construction and operation of a 100 kton liquid Argon TPC certainly represents 
a technological challenge at the present state of knowledge of the technique. For it to 
become a realistic option, we stress the importance and necessity of further R&D and of 
dedicated experimental measurement campaigns. At this stage, we intend to pursue our 
investigations on a ton-scale prototype based on the novel double-phase readout imaging 
method. Options to operate small devices in the J-PARC neutrino beam are being assessed 
in parallel. In addition, we have started to address the possibility of an "intermediate" 
prototype, presumably in the range of 1 kton of mass, based on a similar but scaled down 
design of the potential 100 kton detector. 
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